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Abstract 

The thcrmodynamical properties of a system of two coupled harmonic oscillators in the presence 
of an uniform magnetic field B are investigated. Using an unitary transformation, we show that the 
system can be diagonalized in simple way and then obtain the energy spectrum solutions. These will 
be used to determine the thermodynamical potential in terms of different physical parameters like 
the coupling parameter a. This allows us to give a generalization of already significant published 
work and obtain different results, those could be used to discuss the magnetism of the system. 
Different limiting cases, in terms of a and B, have been discussed. In fact, quantum corrections to 
the Landau diamagntesim and orbital paramagnetism are found. 
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1 Introduction 



Since the pioneering work of Landau in 1930 [1], orbital magnetism of electron gases has been the 
subject of considerable attention, especially during the last decades with the advent of experimen- 
tal opportunities, more precisely with the availability of two-dimensional electronic devices, quantum 
boxes, or mesoscopic finite-size objects. One can find in [2] or [3] a good account of the theoreti- 
cal investigations on the subject, especially from a semiclassical point of view. Fore a more recent 
developments, we cite the book [4] 

Fukuyama group has been developed an amount of papers dealing with different features of two- 
dimensional systems. Among them, we cite the reference [5] where the magnetization of such systems 
in external potentials [5] are studied. In fact, the magnetic field B and temperature T dependence of 
the magnetization is calculated exactly. It is found that the magnetization is well defined in the limit 
of vanishing B as well as in the limit of T = 0, which showing a large fluctuation at low T as B is 
varied. It is shown that this fluctuating magnetization tends to the Landau diamagnetism at higher T 
or by ensemble averaging. Subsequently, other exciting investigation has been reported on [6], which 
concerned the spatial distribution of electric current under B and the resultant orbital magnetism for 
the present system under a harmonic confining potential V(r) = rriLOQr 2 /2 in various regimes of the 
couple (T,B). As an interesting result is that the microscopic conditions for the validity of Landau 
diamagnetism are clarified. 

The system studied in [6] has been considered from another point of view. More precisely, a coher- 
ent states approach is used to investigate its basic features [7]. In fact, the corresponding expressions 
for the thermodynamical potential and magnetic moment are determined. These are exact, in con- 
trast to those in reference [6], and the results yielded a full description of the phase diagram of the 
magnetization. The derivation crucially rests upon the observation that the Fermi-Dirac function is 
a fixed point of the Fourier transform. Exact series expansions ensue by simple application of the 
residue theorem. The related physical quantities are obtained and different discussions are reported in 
terms of the natures of T and B. These concern the thermodynamical potential, the orbital magnetic 
moment, the subsequent magnetic susceptibility and the average number of electrons. 

On the other hand, the problem of two coupled harmonic oscillators living on two dimensions was 
investigated at different occasions where several papers are developed by Kim group, for a short list 
we cite [8]-[15]. Furthermore, the quantum mechanical of such systems on the non-commutative plane 
has been studied as well [16] where different quantum corrections to the original work [8] are obtained 
and their interpretations are given. 

After mentioning the above results, an interesting question arises immediately that concerns other 
features of two coupled harmonic oscillators. Specifically, it is possible to study the thermodynamical 
properties of such systems in the presence of an uniform magnetic field. The answer will be the subject 
of the present paper where interesting results will be derived and discussed. In fact, we will show how 
to use the machinery developed by one of the present authors in the basic reference [7] to analyze the 
magnetization of the system. 

More precisely, we develop a theory that analyzes the basic features of two coupled harmonic 
oscillators under the magnetic field. In doing so, we inspect two already published works [7, 16] to 
generate a full description of the present system from thermodynamical point of view. Actually, this 
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can be done by the help of the energy spectrum solutions. To derive them, we make use of an unitary 
transformation that leads to a solvable Hamiltonian of the system. 

Subsequently, we present two ways to evaluate the thermodynamical potential. Indeed, from the 
Berezin-Lieb inequalities and after determining some physical quantities, we discuss different limiting 
case in terms of the involved physical parameters. These lead to end up with interesting results and 
in particular we show that the average number of electrons behaves like in the inverse of squared 
magnetic field for the infinite coupling limit, i.e. a — ► oo. In this situation, the system behaves like a 
quantum Hall effect one [18]. On the other hand, quantum corrections to the orbital paramagnetism 
and Landau diamagnetism are obtained. More importantly, we notice that by switching off a in our 
analysis, we recover already published work [7]. 

Furthermore, we give an exact formula of the thermodynamical potential with the help of some 
well-known relations and discuss different issues. Indeed, applying the Fermi-Dirac trace formulas, 
we explicitly derive the average number of electrons and the magnetic moment. Finally, we how that 
they can be reduced to the standard expressions by taking into account the the liming case a = 0. 

The present paper is organized as follows. In section 2, we formulate our problem by establishing 
the necessary materials to deal with our task. In section 3, after making use of an unitary transfor- 
mation, we introduce an algebraic method to derive the energy spectrum solutions. Their underlying 
properties will be discussed by considering four limiting cases. We construct the coherent states for the 
present solutions and show that they are coupling parameter dependent in section 4. These will serve 
as tools to determine explicitly different physical quantities and in particular the thermodynamical 
potential in section 5. It will be obtained by adopting the Berezin-Lieb inequalities in the first stage. 
This allows us to give different discussion and end up with interesting conclusions. However, in section 
6, we consider another approach based on the Fermi-Dirac trace formulas to give the exact form of 
the thermodynamical potential. Finally, we conclude and give different perspectives. 

2 Formulating the problem 

We start by formulating our problem by setting the needed tools for doing our task. This can be 
done by establishing a mathematical formalism governed by a Hamiltonian describing a system of two 
coupled harmonic oscillators in two dimensions. Subsequently, we submit the system to a constant 
magnetic field and analyze its behavior. In doing so, we determine the energy spectrum through an 
algebraic method after making use of an unitary transformation. 

2.1 Coupled harmonic oscillators 

We consider a system of two coupled harmonic oscillators of mass (mi, 772-2) and living on the plane 
(Xi, X%). This can be described by a Hamiltonian as sum of free and interacting parts, such as 



where C\, C2 and C3 are three constant parameters. Note in passing that, the involved parameter 
can be fixed according to the nature of the system. On the other hand, (1) has been investigated 
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for different purposes, for instance one may see reference [8], which has been generalized to the non- 
commutative geometry case [16]. 

As claimed before, we are wondering to study the magnetization of two coupled harmonic oscillators 
under an uniform magnetic field. To achieve this goal, we generalize the system governed by (1) to 
another one of Hamiltonian 

Pi = — + — + - (CiX? + C 2 Xl + C 3 X X X 2 ) (2) 
2roi 2ma 2 v 1 2 ' 

where IIi and II2 are the conjugate momentum. They can be simplified by choosing an appropriate 

gauge. Indeed, in the symmetric gauge 

A = ^(-X 2 ,X l ) (3) 

they are given by 

eB eB 
n 1 = p 1 -—x 2 , U 2 = P 2 -—X 1 . (4) 

Using these to map (2) into the form 

Pi = P- + p- + \ {D x Xf + D 2 Xj + £3*1*2) + \ (oj 1 P 2 X 1 - u; 2 PiX 2 ) (5) 
Zvtl\ Zm 2 Z Z 

where the new constants D±, D 2 and D 3 read as 

D 1 (B) = C 1 +rrw4, D 2 (B) = C 2 + m x J{ , D 3 = C 3 (6) 
with the cyclotron frequencies 

eB eB 

ui c = , u 2c = . (7) 

m\c m 2 c 

Clearly, by comparing and forgetting about different involved constants, we notice that the third term 
makes difference between (1) and (5). This in fact will play a crucial role in the forthcoming analysis 
and allow us to derive different results. 

It is convenient to introduce new phase space variables, which can be done by rescaling those 
appearing in (5). Indeed, one can define the positions as 

*i=(^y*i, x ,={^fx 2 ( 8) 

\m 2 J \mij 



which obviously lead to the momenta 



1 1 
m 2 \ 4 r> ( m\ \ 4 



Pi = — Pi, P2 = — P 2 - (9) 
\mij \m 2 ' 



Replacing all, we show that (5) becomes 

H 2 = — (pf + p\) + - (dix\ + d 2 x\ + d 3 xix 2 ) + ^ (xip 2 - x 2 pi) (10) 



where we have set different constants as 

1 1 

d\{B) = D\ J — J 2 d 2 (B) = D 2 (^) , d 3 = D 3 (11) 

with unique mass m = (mim 2 )2 and the cyclotron frequency uj c = {w\uj 2 )^ = Consequently, (10) 
is showing up an extra term, which is nothing but the angular momenta and the first is similar to (1). 
Therefore, it will be of interest to deal with such system and underline its physical properties. 
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2.2 Unitary transformation 

According to the expression form (10), it appears that getting the energy spectrum solutions is not 
a easy task. However, we can overcome such difficulties by adopting an appropriate approach. More 
precisely, we proceed by making use of an unitary transformation, such that new phase space variables 
can be defined by 

y a = M ab x bl p a = M ab p b (12) 

where the matrix M ab 

(cos | — sin I \ . . 

■8 ) ( 13 ) 

sin g cos 2 / 

is an unitary rotation with the mixing angle 6. Inserting the mapping (12) into (10), one realizes that 
should satisfy the condition 

tang = - d \ (14) 
d 2 — d\ 



to end up with a factorizing Hamiltonian. It is 

H3 = i + ^ + \ ^ + e ~ 2ay ^ + f {V1P ~ 2 ~ VW ~ l) (15) 
where k and a are given by 



fc = y^ 2 -f, ea = d 1 + d 2 + ^(d,-d 2 y + di (l6) 

and the condition Adid 2 > d^ must be fulfilled. Note that, H% has a form similar to two-dimensional 
Landau Hamiltonian in the symmetric gauge. Obviously, they coincide in the case of without coupling, 
namely a = 0. 

Before proceeding further, we conclude by citing some interesting remarks. In doing so, let us 
return to H% and define two operators as 

H ° = ~L, ^ + ^ + \ (g2 °^ + e ~ 2ay ^ ' Ls = {yw ~ 2 ~ mpi) (17) 

where Hq also can be separated into two commuting parts 
Firstly, one can see that the decoupled Hamiltonian 

can be recovered by taking a = 0, which corresponds to the solution d\ = d 2 and d% = 0. Secondly, it 
is interesting to note that (19) can be derived by a canonical transformation only from 

H = H 1 +H 2 (20) 

as it is pointed out in [8] and subsequently in [16]. 

According to the above statements, we can rearrange H$ in an appropriate form. This is 

H 3 = e~ a Hi + e a H 2 + L 3 (21) 

which will be used to tackle different issues in the forthcoming analysis and in particular the magne- 
tization of the present system. This mapping will be helpful in sense that the corresponding energy 
spectrum solutions can easily be obtained as we will see soon. 
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3 Energy spectrum 



As far as the eigenvalues and eigenstates are concerned, we adopt an algebraic method based on 
different operators in terms of the phase space ones. This will allow us to obtain the solutions and 
investigate their underlying properties. 

3.1 Algebraic analysis 

It is clear that Ti is a Hamiltonian of two decoupled harmonic oscillators. Thus it can simply be 
diagonalized by defining a set of creation and annihilation operators. They are given by 



k c* 2 a_ _ -j- / k c* % 



e2?/H e 2tL a J = W^^ e2 y- e 2 a (22) 

where the new frequency is 

, v (U^-dW 1 l~k , s 

- (B) = Hi^) = V™ (23) 

They satisfy the usual commutation relations 

>',.«;: »,j (24) 

and obviously other commutators vanish. It is easy to show that Ti can be mapped in terms of ai and 
a] as 



H = hw {a\a x + a\a 2 + l) . (25) 



According to (25), it is not hard to derive the corresponding energy spectrum solutions. This can 
be done by solving the eigenvalue equation 

H\n 1 ,n2, a) = £ nim |rai, n 2 , a) (26) 

to get the corresponding states 

K,n 2 ,a) = ^4^-|0, 0, a) (27) 

as well as the energy spectrum 

£ ni ,n 2 = ( n i + n 2 + 1) • (28) 

Due the fact that there are mappings between different Hamiltonian's, one can build other solu- 
tions. In particular, the spectrum of Hq can easily be deduced from above as 



(29) 



To get that for H3, we need to diagonalize the angular momentum. In doing so, we define two sets of 
operators where the first one is 

0,9 = 71 + ia2 ^ ' 0,9 = 71 ^ ~ ia ^ ^ 
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and the second reads as 



a>d = "i - ia2 ) ' a d = -Jk ( a i + ia £) ■ ( 31 ) 



They are showing 



cig, o)g 



a d , 4 



1 (32) 



and different commutation relations are nulls. One can note that there is a conservation of number 
operators, such as 

N 1 + N 2 = N d + N g . (33) 

where we have N\ = a\a\, N 2 = a\a2, N d = a d a d and N g = a g a g . Now we express the phase space 
variables in terms of the new operators to end up with a quantized angular momenta. This is 

L 3 = 2h(N d -N g ). (34) 

Obviously, its eigenvalues are 2h (n d — n g ) and the corresponding eigenvalues are forming a common 
basis of L 3 and TC. 

Finally, we settled all ingredients to derive the energy spectrum solutions of H 3 . Indeed, starting 
from the above results, one can see that (15) becomes 

H 3 = (hwe a + Huj c ) N d + (hwe~ a - hw c ) N g + hw cosh a. (35) 

To write #3 in compact form, it is convenient to introduce two new frequencies in terms of the former 
ones. These are defined by 

u>+(B, a) = we" + u c , U-(B, a) = uoe~ a — u c . (36) 

They are showing a strong dependence to a and therefore generalize the standard results [7]. Now 
returning to map H 3 as 

H 3 = h{uj + N d + oJ^Ng + w cosh a) . (37) 

Solving the eigenvalue equation, we can easily derive the energy spectrum solutions. Thus, the 
eigenvalues take the form 

Ez,n d ,n g = ft {u+n d + uj^rig + u cosh a) , n g ,n d = 0, 1, 2, • ■ ■ (38) 

and the eigenstates are given by 



(4W 



\n d , n g , a) = — 10, 0, a). (39) 

\Jn d \n g \ 

It is clear that the results obtained so far are a-dependent. This in fact makes difference with respect 
to the standard results obtained by analyzing the Fock-Darwin Hamiltonian [7], which obviously can 
be recovered by setting a = 0. At this stage, one may ask about the relevance of such coupling 
parameter and the answer will be given in the forthcoming sections where interesting results will be 
derived and different discussions will be given. 
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3.2 Underlying properties 

In investigating the underlying symmetry of the system, one can study the properties of quantum 
numbers pairs (n^,n 9 ). However, these may not provide simple hints on the ordering of the energy 
E a ,n d ,n g with the exception of four limiting cases related to the nature of the coupling parameter and 
the magnetic field. 

3.2.1 Weak coupling case 

To characterize the system behavior, we consider the first case that corresponds to the limit a — > 0, 
which means that the coupling is not strong enough between two oscillators. This is the case for some 
physical phenomena. Therefore, we can make different approximations to approach our findings to 
well-know and significant results. 

By taking the limit a — > and after a simple calculation, we show that the energy spectrum can 
be approximated by 

E3,n d ,n g \a^0 ~ K [u(n d + n g) + ( auJ + u c ){n d - n g ) + u] . (40) 
We can bring this to an appropriate form by defining new quantum numbers. They are 

A=-^-£, ^ = ^V^- ( 41 ) 

Thus, one can rearrange (40) as 

UJ 



Ea,\,£,\a->0 ~ 2f>< 



uj\ + (auj + uj c )i + - . (42) 



It can be identified to the eigenvalues of the Fock-Darwin Hamiltonian, which can be obtained from 
(38) by taking a = 0. They are 

E3,n d ,n g | a=0 = ^ (^+n d + UJ-U g + u) . (43) 

Now observing that the following correspondence (cj+,a;_) — > (w,aw + u c ). This tells us that (42) 
can be used to analyze the thermodynamical properties in similar way to that has been done in [7]. 
Moreover, it shows how one can generalize the Fock-Darwin Hamiltonian to another one where the 
interaction still surviving. 

At this stage, we can further discuss (42) by inspecting two other limits in terms of the field. In 
doing so, we suppose that the cyclotron frequency is much smaller than the frequency uj, i.e. uj c <C oj, 
thus we have 

« 2Huj [X + a£ + . (44) 

According to this, two conclusion can be deduced here. Indeed, firstly we still have a generalized 
Fock-Darwin Hamiltonian but its frequencies are changed now to (w,qw). Secondly, without coupling 
we recovers one-dimensional harmonic oscillator of eigenvalues 

E x \ a=0tB _ t0 -2huj(x+^\ (45) 



whose frequency is u = \ ^r, which means that we are in the conditions d\ = d,2 and = 0. 
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Now, let us treat the second consideration that is the strong magnetic filed case. In fact, this 
equivalent to uj c 3> ui and leads 

#ela-+o,£-»oo, ~ 2hw c £. (46) 

It can be interpreted as the squared energy spectrum of the massless Dirac fermions in graphene under 
an uniform magnetic field. Fore more detail, we cite for instance [17] and reference therein. 

3.2.2 Strong coupling case 

It is immediate and natural to ask about what happens if the coupling is strong enough and the 
corresponding limit cases of magnetic field. 

The above inquiry can be answered by examining the limit a — ► oo. Thus, returning to (38) to 
show the result 

Ea,n d ,n g \a— >oo 

Again this can be approximated further by taking other limits. Indeed, focusing on the case lj c <C oj 
or ijj c <C uie a , it is straightforward to obtain 

E nd ~hio(n d + e a (48) 

which is the energy spectrum of harmonic oscillator of frequency u;e a . However for oj c S> uj, there is 
nothing to say and therefore (47) remains as it is because we can not make comparison. 

In summary, according to the above results we conclude that the coupling parameter a is interesting 
parameter of the present theory. In fact, it can be adjusted to recover different models those used to 
deal with different issues in physics. 



we 



n d + - ) + lo c (n d 



11, 



(47) 



4 Realizing the coherent states 

The forthcoming analysis requires a powerful tools. More precisely, one way to determine the ther- 
modynamical potential is to use the coherent states approach. Thus, for the neediness, we follow the 
standard method to realize them in terms our language and show their dependence to the coupling 
parameter. In fact, we will use the same steps traced in [7]. 

The fact that the eigenstates issued from the algebraic method are just tensor products of Fock har- 
monic oscillator eigenstates allows one to easily construct the corresponding coherent states. Indeed, 
in a standard way, we have 

| z d , z g , a) =| z d ) <8> | z g ) = exp 

In terms of the creations operators, we have 
| z d ,z g ,a) = exp 

The above normalized states, should obey some of the usual properties. Indeed, it is easy to verify 
the eigenvector property, such as 

a d , \ z d ,z g ,a) = z d \ z d ,z g ,a), a g \ z d , z g ,a) = z g \ z d , z g ,a). (50) 



-U\z d \ 2 + \z 9 f) 




rid, rig. 



(I* 



2 + N 2 ) 



z d a\+z g a 



•)"g 



0,0, a) 



(49) 
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As far as the action identity is concerned, one can obtain the relation 



H 3 (z d ,z g ,a) = (z d ,z g ,a \ H 3 \ z d ,z g ,a) = h(iv + \z d \ 2 + LO_\z g \ 2 + to cosh a) (51) 

where the function H 3 (z d , z g , a) is called lower symbol of the operator #3. It will plays an important 
role in the present context. The resolution of the identity reads as 

l= — I I z d ,Zg,a)(z d ,z g ,a \ d 2 z d d 2 z g . (52) 

where the last property is also crucial in our context. 

For any observable A with suitable operator properties (traceclass, • • • ), there exists a unique 
upper (or covariant) symbol A{z d ,z g ) defined by 

A=-y I A(z d ,z g ,a) I z d ,Zg,a)(z d ,z g ,a \ d 2 z d d 2 z g . (53) 

As a straightforward illustration, we consider the upper symbols for the number operators. Hence, 
one can show 

N d (z d ,z g ,a) = \z d \ 2 -1, Ng(z d ,z g ,a) = \z g \ 2 - I. (54) 
Clearly, the upper symbol for our Hamiltonian (37) takes the form 

Hs(z d , z g , a) = h (tu + \z d \ 2 + ui-\z g \ 2 — uj cosh a) . (55) 

To setup all what we need for our task, we recall an useful trace identity for a given traceclass 
observable A. This is 

TrA = -=[ A(z d ,Zg,a)d 2 z d d 2 z g = [ A(z d , z g , a) d 2 z d d 2 z g (56) 

7T JC 2 7T JC 2 

where the symbol function A is 

A(z d , Zg, a) = {z d , z g ,a\A \ z d , z g , a). (57) 

We close this part by noting that all involved quantities are a-dependent. Obviously, the standard 
results can be recovered by switching off the coupling parameter [7]. On the other hand, we will see 
how the above materials can be employed to deal with different issues and in particular determine 
the thermodynamical potential. This will be done by adopting two methods, which concern the 
Berezin-Lieb inequalities and Fermi-Dirac trace formulas. 



5 Berezin— Lieb inequalities 

Having derived and settled all necessary tools, we now show that how they can be used to study the 
magnetism of the system under consideration. In doing so, we start by defining the physical quantities 
those will be discussed in the present context. One way to do so is to evaluate the thermodynamical 
potential, which can be done, in the first stage, by adopting the Berezin-Lieb inequalities. Subse- 
quently, we treat the asymptotic behavior of the obtained results by considering the liming cases of 
the coupling parameter. 



9 



5.1 Physical quantities 

The magnetism of the model under hand can be investigated by adopting the standard method of 
statistical mechanics. This will be done by making use of different approximations to simplify our 
problem. In fact, we begin by assuming that the total number (N e ) of electrons is large enough for 
making no appreciable difference between a grand canonical ensemble and a canonical one. 

On the light of the above considerations and obtained results, we proceed by using the magnetic 
moment M definition. This is 

where thermodynamical potential f2 can be obtained from the partition function. In terms of our 
model, it is 

n = -^Trlog [1 + exp{-/3(tf 3 - p)}] (59) 
as usual we have set (3 = 1/(^7). According to (58) and (59), we show the result 

M = -2fi B Tr ( Nd ~¥ • (60) 

Replacing H 3 by its expression and tracing to end up with 

M = -2/x B y (61) 

n d ^=o 1 + K ± exp{^(w + n d + a;_n s )} 

where /i_B = He/(2mc) is the Bohr magneton and k± are given by 

k±(B,T, a) = exp \fl(fJ> db tkJ cosh a)]. (62) 

On the other hand, the average number of electrons can be evaluated by introducing the Fermi 
distribution function. That is 



Therefore, in our case we have 



(N e )= f(E 3 , nd ,n 9 ) = Tvf(H 3 ) = -d fi n (64) 

n d ,n g =0 

which is showing that there are two possibilities to get (N e ) either summing all distributions (63) or 
deriving Q with respect to chemical potential. Clearly, to go further in evaluating different physical 
quantities, one should explicitly determine f2. 



5.2 Calculating the thermodynamical potential 

The thermodynamical potential is very much needed to describe the quasi-classical behavior of present 
system. This can be calculated by adopting some technical methods like for instance the Berezin-Lieb 
inequalities. In fact, it is based on some general statement that is for any convex function g(A) of the 
observable A, one can write the inequalities 

\! g{A)d 2 z d d 2 z g <Tig{A)<\[ g(A)d 2 z d d 2 z g (65) 
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where the lower and upper symbol functions (A, A) are defined before. This tells us that knowing 
inferior and superior boundaries of a given observable, one can derive its trace. 

At this level, we have all ingredients needed to do our task. Indeed, an straightforward application 
of (65) gives the result 

- 7T2 / lo § f 1 + ex P{-/W - »)}] d 2 z d d 2 z g < Q 
p^ J : ••• L J 

^ "A f lo § t 1 + ex P{-/W - /*)}] d 2 z d d 2 z g . (66) 
After mapping (51) and (55) into (66), we end up with 

^ f'OO roc 

— — j dud / du g log [1 + exp{— f3(huj + Ud + hw-u g — tuv cosh a — //)}] < Q 
P Jo Jo 

roc roo 

< — — / (in^ / du g log [1 + exp {— (3(ftw + Ud + ftw^u g + huj cosh a — fi)}] (67) 
P Jo Jo 

where we have set Ud = \zd\ 2 and u g = \z g \ 2 . The solution can be obtained by making some rearrange- 
ment followed by an integration. Indeed, by changing variables as 

u = f3h(u + Ud + U-Ug), v = f3hu> + Ud (68) 

we show that, in terms of the parameters k± (62), (67) becomes 

$(/e + ) < Q < $(«_) (69) 

where the function $ is given by 

= — — — / du. (70) 

Actually, the problem of determining Vt is restricted to find the solutions of such integral. This can 
be done by defining a new parameter 

A = uj + u)^ = uj 2 — uj 2 — 2ujuj c sinh a (71) 

and distinguishing between the sign of k± to end up with the solutions 

*f^_J mw^i-K), k<i 
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K > 1 



where we have introduced here the function F s of the Riemann-Fermi-Dirac type. For a given variable 
z, it reads as 



oo 

Z 



m 

m=l 



Note in passing that $ is depending to the sign of A as well. Discussions about such matters will be 
reported next. On the other hand, the above results can be discussed by separately considering the 
high and low temperature regimes. 
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5.2.1 High temperature regime 

Having the expression (69) together with (72), one can introduce an appropriate approximation to 
further simplify the form of O and derive interesting results. To achieve this goal, we can analyze two 
liming cases of the temperature of the present system. 

We start our analysis by considering the high temperature regime that corresponds to the condition 
\fx ± 7ia;cosha| <C k B T. By taking into account, we find 



n 



^bT ( k B T\ 



A 



Fs(-1) 



-0.901543 



knT 



(k B T\ 2 



x \ h J 



(74) 



This is a nice form that can be further discussed. Recall that the involved parameter A is magnetic 
field and a-dependent, which gives a generalization to the already obtained in [7]. Obviously without 
coupling, they coincide. 

On the other hand, one can report different discussions related to the above form of SI in terms of 
the coupling parametera. With these we can show what makes difference with respect to the standard 
case, i.e. c\ = C2 and C3 = 0. By doing this, we can summarize the following results: 

• By inspecting the form of f2, one can immediately notice the first general result. Indeed, by 
considering a negative A we end up with a positive £1, which can not be obtained from the 
standard results [7]. 

• (74) is magnetic field dependent as well and therefore the present case exhibits an magnetism 
behavior. This statement can be confirmed by explicitly determining the magnetic moment and 
susceptibility. 

• The easiest way to obtain the magnetic moment is that one can require for instance the following 
configuration: 

u — > lo c , sinh a = finite number. (75) 



After giving quick conclusions by looking at the form (74), now let us be much more accurate and 
derive explicit results. Indeed, after a straightforward calculation, we show that the magnetic moment 
takes the form 



fk B Te\ 
M = 1.803086 k B T — 

\ Xnmc J 



B 

2^ 



C2 + m\U)\ _|_ c l + \ 



mi 



m 2 



sinh a 



ujmc 2Blo c ( C2 + 



This allows us to end up with the susceptibility 



X = -3.62172 k B T 



fk B Te\ 
1 kmc J 



8m 



4cic 2 



mi 



sinh a + 



+ 



m 2 



4cic 2 



\mi rri2 



+ 1 



(76) 



(77) 



It is clear that x is behaving as a linear function in terms of the hyperbolic function sinha. This 
results in fact is showing the difference with respect to the case without coupling where there is no 
susceptibility and therefore no effect is obtained at high temperature regime. 
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5.2.2 Low temperature regime 



To accomplish our analysis in terms of temperature, we discuss the last case. This can be achieved 
by considering the more realistic case, which is \i ^> Hujsmha and 11 S> k B T . With these, we will be 
able to derive interesting results and deduce different conclusions. 

After considering the above two limiting cases, we show that (72) can be written as combination 
of three parts. This is 



%) = i T - + 5± 



(78) 



where different terms are given by 



A(B,T,a) 



j(B,T,a) 



S±(B,T,a) 



-A 

2A 3 \h 

Hui cosh a 
2A 

kuT f knT 



A 



M 2 , 1 2 ,2 ^ fk B T\ 

v + r cosh a+ Y{—) 

2 

Fz(— exp [— ± Huj cosh a)]). 



(79) 



n 



This will 



According to these functions, we notice that f2 is in the interval [A + S + — — , A + S- + 
be used to derive different results in the present context. 

We can go further by making an important assumption. In fact, we restrict ourselves to the 
condition 

In this situation, one can see that S± is reduced to 

2 



S± w S (B,T,a) 



k B T ( k B T 



A 



Hurt 



Moreover, taking into account the above limiting cases, we show 

A 



\A + S \ 



0. 



Combining all to end up with the form 

2 



ix 1 



2 I 3 V Hu, 



Jo 



2 V.2 , * ( k B T 

uj cosh a H — — 

3 V Hujq 



Wig)* 



(81) 



(82) 



(83) 



This in fact can be used to deduce different physical quantities. In particular, we evaluate the average 
number of electrons to obtain 



(jv e >(B,r,a)«i(^ ! 



1 



77+ — 



Huj\ ' 



V J 



vr 2 (k n T 



+ 



k B T 
A* 



F-> 



(84) 



As far as the magnetic moment is concerned, one can obtained a complicated form. This is due to 
the fact A and uj are magnetic field dependents. However, we can get more information by inspecting 
some limiting cases. These will also offer for us a way to emphasis what makes difference with respect 
to other approaches and in particular [7]. This can be done by discussing the nature of the coupling 
parameter involved in the game. 
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5.3 Asymptotic behavior 

Having derived a general expression of the thermodynamical potential, one can ask about further 
simplifications to characterize the system behavior in some special cases. More precisely, how the 
above results can be approximated by inspecting the limits: a = 0, a <C 1 and a — > oo. The reply of 
such question is the subject of the next investigations. 



5.3.1 Without coupling 

It is natural to ask about the case a = 0. To reply this inquiry, one can return to the former analysis 
to show that the different quantities given in (79) can be restricted to the functions 



A(B, T, 0) 



2 



JJ_ 



+ 



u_\ 2 | 7T 2 /k B TY 



A, . Huj 

j(B,T,0) = — 



(/^o) + 3 (w ) + 6 



k B T 

fkOn 



(85) 



S ± (B,T,0) 



k B T [ ^) F 3 (- exp [-p(ji ± fiw)}). 



They show that the thermodynamical potential lies in [A + 5+ — —, A + S+ + ^1 . Note that, these 
exactly coincide with those obtained by analyzing a confined two-dimensional system in the presence 
of an uniform magnetic field [7] . 

To reproduce most of results derived in [7], one can inspect (85) by making an approximation. 
That is S± can be replaced by 



S = k B T 



( k B T 



Fa 



to end up with the form 



2 } 3 \hw J 



U \ 7T 
UJ ) 3 



/ k B T\ 2 ) +k T (k B T\ 2 F 



(86) 



(87) 



Therefore, the average number of electrons is given by 

2 / \ 2 



-Pp. 



as well as the magnetic moment 



The corresponding susceptibility Xp rea d as 



(89) 



(90) 



One can also inspect other approximations. Indeed, by requiring that /i 3> k B T and fi <C fkv, we 
show that (88) can be written 

2 \ 

The above derivation show that our results are general in sense that after making appropriate choices 
one can recover already significant published works. 



{N e 



(91) 
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5.3.2 Week coupling 



We start our analysis by dealing with the first case that corresponds to «< 1. Clearly, an expansion 
of different quantities entering in the game is very much needed. Indeed, by taking the first order of 
a, we can approximate (83) as 



n 




\oj J 3 \frw J 



+ k B T 



\hw J 3 



-/V 



(92) 



where we have set uJq = oj 2 — uj 2 . 

Obtaining (83), it is worthwhile to ask about the related physical quantities to characterize their 
behaviors in terms the coupling parameter for the present case. Using the former definitions to show 
that the average number of electrons is 



(N e )(B,T,a) 



r + — 



v J 



vr 2 (k B T\ 2 fk B T 



F, 



1 + 



2uJUJ r 



-a 



(93) 



On the light of the assumptions fi S> k B T and \jl 3> hiv, we obtain 



{N e )(B,a) 



2 \hw J 



+ lOLO c 



a. 



(94) 



Clearly, the second term in right hand is appearing a correction to the average number of electrons. 
This is agreed by canceling the coupling to recover the standard result (91). 

Now let us investigate the magnetism in such case. Indeed, a straightforward calculation gives the 
magnetic moment as 

,2- 



M(B,T, a) 



Huiq 



flUJ c 1 + 



2lulu c 



-a 



W + L0 r 



n(B,T, 0) 



(95) 



where Q(B, T, 0) is the thermodynamical potential corresponding to the standard case, i.e. a = 0. 
Considering /i k B T and [i S> hu, we show 



3 /jcJq \ m^o / 



(96) 



One important thing should be noted here is that the magnetic moment is behaving like a linear 
function in terms of a. Obviously, without coupling we end up a null magnetization. This means that, 
we have like phase transition from coupling to decoupling system. This point might be investigated 
further to deal with other issues in statistical physics. Furthermore, by carefully identifying (89) to 
(96), one can fix a to reproduce the orbital paramagnetism. Indeed, the solution can be written as 



a 



2hw\\i B ijj 
UJ 2 + u> 2 



B. 



(97) 



On the other hand, one can also make another choice of the coupling parameter to get interesting 
result (89). In fact, here also one can reproduce the Landau diamagnetism. 
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5.3.3 Strong coupling 

To complete our analysis we consider the last case that is the strong coupling limit. This of course 
will shine light on the system behavior at such case and therefore allow us to get more interesting 
results. To clarify this, we take the limit a — ► oo to obtain 

A- Ml-)', £-£[T>V S ± =0 (98) 



U) C J 2 12 \ui, 

which leads to the thermodynamical potential 



to, 



2 



It is clear that the average number reads as 

which behaves as the inverse of magnetic field. It seems that (100) is sharing some common features 
with the quantum Hall effect results [18]. Indeed, one has to recall that the filling factor is defined as 
the ration between {N e ) and the quantized flux. More precisely, we can write 

¥• "*=¥■ (101 » 

where S is the system area and (j>Q = —. Clearly, we can adjust all parameter to show that effectively 
we have something related to the quantum Hall effect. 

In summary, the Berezin-Lieb inequalities are a powerful tools one can use to study the thermo- 
dynamical behavior for a given system. As we have seen so far, a straightforward application of such 
approach allows us to derive different interesting results. On the other hand, as we claimed before 
there is another way to do so and this will be tackled next. 



6 Fermi— Dirac trace formulas 

As we claimed before, we use the second method to explicitly determine the exact expressions for 
the thermodynamical potential. This is in fact based on the Fermi-Dirac trace formulas, which does 
not include include approximations in the derivation of £1 and therefore makes difference with respect 
to the Berezin-Lieb inequalities. Subsequently, we restrict ourselves to the evaluation of the average 
number of electrons and magnetic moment as well as their expressions at zero coupling. 

6.1 Exact expressions of f) 

Using the machinery developed in the reference [7], we can derive an exact form of 17. This can be 
done by making an straightforward application of the Fermi-Dirac trace formulas, in particular (139) 
and (141) in the appendix. Using (140) to define a function Q(k) in terms of language as 

9(fc) = Tr [ e -(i*+i)f *&1 (102) 
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where the Hamiltonian H3 is given in (15). After replacing #3, we end up with 

e(fc) = Tr 



-(ik+l) £§■ {ui + N d +ui- N g +u) cosh a) 



which can be written as 

Q(J^ — g—(*fc+l)f Neosho 



1 



1 



1 _ g -(ifc+l)f 1 _ e -(ifc+l)§fiw- 



(103) 



(104) 



Now let tackle our problem by writing the Fourier integral representation for the thermodynamical 
potential. This is 

r+00 g— (ik+l) & (fw cosh a— fi) 



ft 



2 cosh 5 k 



1 



ik + l J \ i — e -(ifc+i)f / V 1 



-(ife+l)|na>- 



dfc. (105) 



This integral is given as a series by using the residue theorem. One can easily see that the numbers 
(2m + m G Z are simple pole of cosh and i + y^— y < 
poles of Q(k). Now we can consider two case, the first one where 



4™ ^ - + 4p» m G ^* are simple or double 



a G 



-00, log 



fvjj 



— — \ / ( —\ 



\hwJ 



1 



u 



log 



fuj V \hw 



1 L+OO 



(106) 



and here we take an integration path lying in the lower half-plane and involving only the simple poles 
(2m + m < 0. It leads to the result 



Q{B,T, a) 
In the second case where 



. 00 
4/3 2 



("1) ? 



m sinh (f hw + m) sinh (^HuJ^m) 



(107) 



a G 



log 



\ftjbj' 



1 Jos 



— + A / ( — \ 



\fVjjJ 



1 



(108) 



an integration path in the upper half-plane is chosen. It encircles all the other poles: (2m+l)z, m > 0, 
i + Airm/ ((3fujo + ), i 4- Anm / (j3ftw _) , m G Z*. We present the result in a manner which will render 
apparent the various regimes 



ft 



(ft L + ft i) + 



02 



+ 



2vri| 



o_i(j) 4- ^ o_i [(2m + l)i] + ^ 

m>l rn±=£0 



4-7T 



(109) 



where Q,i(B,a) is given by 



and Oqi(-B,T, a) reads as 



ft 7 



24A 24 \w 2 -w c 2 - 2wu^ sinh a 



(110) 



n 



01 



6 



%\J uj 2 — ujI — 2ujuj c sinh a 



4-^ 



fty — u> 2 — 2ujuj c sinh a 



(111) 
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Now let us consider an approximation such that 
expansion 

1 



2lulu c 



sinh a <C 1. This allows us to write the 



n — 2uJUJ r sinh a 



~ 1- 







2lULU, 



sinh a 



in the first order. Consequently, we find This gives Ql(B,o) as 



n L (B,a) 



1L. 

24 



u 



+ 



12 



coujr sinh a. 



(112) 



(113) 



This result can interpreted in different ways. Indeed, if we forget about the a appearing, we can reach 
the same conclusion as in [7]. Indeed, the first term is at the origin of the Landau diamagnetism and 
gives the susceptibility 



1 



XL 



MB V 



1 



(114) 



where the coefficient Dn 



can be interpreted as the density of states at Fermi energy. Note that, 



the value of xl is equal to one third of the one Xp found in (89). On the other hand, a can be adjusted 
to get another contribution to the susceptibility. More precisely, we can define a in terms of the inverse 
of squared magnetic field to absorb the term lolo c and therefore get a a correction to the standard 
Landau diamagnetism. Otherwise, we can even reproduce Xp simply by making an appropriate choice 
of a. Indeed, fixing sinh(a) = to end up with (89), where 7 is constant that can be fixed easily. 
This show how the obtained results generals and allow to deduce interesting properties. 
According to (112), fioi(a) becomes 



fioi(a) 



6 



f2o2 (a) reads as 



^02 (a) 



k B T 



3 



1 £ (-1) 



+ 7T 



LULU, 



sinh a. (115) 



4/3 



m=l 



m 



sinh 



2k B T 



m 



sinh 



2k B T 



-m 



(116) 



Note that, when a = we get A = luq, which leads to recover the result obtained in [7]. It becomes 
negligible at low temperature regime k B T <C \l. The sum of Ql and O01 is analogue to the term A in 
(78) and f2o2 corresponds to S±. The last term is responsible for the oscillatory behavior. If lu+/lu- 
is irrational values, we have 



^osc(a) 



1 



-Y 

2(3 



(-iy 



m=l 



111 



8111 



8111 



2/1 

h(uie~ a —u>c) 



irm 



+ 



sin 

2// 



Lue a +u) c 
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sinh 



2k B T 



7rm 



sm 



-TTIJI 



sinh 



2k B T 
h(uje a +LU c ) 



Tr 2 m 



= ^osc(«) + ^sc(a)- 



(117) 



6.2 Average number of electrons 

In this section, we will exploit the formula's (107)-(117) to obtain the exact expressions of the average 
number of electrons and the magnetic moment. We will restrict ourselves to the more realistic case: 
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/i < hw/2. The average number of electrons is easily derived by taking the derivative of — £1 with 
respect to \i. It is found to be 



(N e ) = -d^n L (a) - df,n 01 (a) - d^n osc (a) - M fi+ c (a). 
With the straightforward calculation, we find 

(N e (a)) = (N e ) L + (N e ) 1 + (N e ) 02 + (N e )- c + (N e >+ sc 

where different portions are give by 



(Ne)L 
(N ) 1 
(N ) 02 

(Ne >osc 
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2/1 -rm 
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n{uie a — ujc) 



in=l 



k B T cos 
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(118) 
(119) 

(120) 
(121) 
(122) 

(123) 

(124) 



The above results are general in sense that the standard solutions can be recovered. Indeed, 
requiring that a = 0, we show 
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24 V^o 

1 00 
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fi<i»-)- 



m) 



6.3 Magnetic moment 



The magnetic moment is found by the conversely of the derivative of the thermodynamical potential 
by the magnetic field. Then, we have 



M(B, T, a) = -2(i B 
It is equivalent to 



V dhw c 



, dn 01 (a) \ f dn 02 (a) 



2ji B 



dfkOr 



M{B,T,a) = 2fi B [Ml (a) + M 01 (a) + M 02 (a) + M~ c (a) + M osc (a)" 



(125) 



19 



where different quantities read as 

M L (B,T,a) = -r^-(^ 



hw 2 , 
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+ 7T U? 
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In the case where a = 0, we end up with 
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and, for the irrational case 0J+/0J- Q, 
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(133) 



(134) 



We will not give the expressions of M^ sc in the rational case because the magnetization is a continuous 
function of ui c and its behavior can be fully understood from the irrational one. 

In the end note that, the temperature scale is compared to the two natural modes w± of the system 
and draws three possible intrinsic regimes: high temperature regime k B T > hw+, low temperature 
regime k B T < hw-, and intermediate temperature regime hw- < k B T < hw + . Remember that we 
work in the large electron number region: fi > hw/2. 
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7 Conclusion 



We started by formulating our problem in two-dimensional space where two coupled harmonic oscilla- 
tors living on. Subsequently, we introduced a minimal coupling to generate another interacting system 
that is studied. After rescaling different variables, we showed that it is possible to get a diagonalized 
Hamiltonian. In fact, this is done by making use of an unitary transformation. It was helpful in 
sense that the eigenvalues and their wavefunctions are obtained in simple way in terms of the coupling 
parameter a. 

The fact that the energy spectrum solutions are (a, B) dependent, we discussed their underly- 
ing properties. More precisely, four limiting cases have been investigated, which are week and strong 
parameters (a, B). In particular, we noticed that by fixing a, some model can be recovered, these con- 
cern for instance the Landau Hamiltonian in two dimensions and harmonic oscillator in one-dimension. 
These allowed us to conclude that by adjusting the coupling parameter, one can derive other interesting 
other solutions. 

To investigate different issues related to the considered system, we constructed the corresponding 
coherent states, which are obtained to be coupling parameter dependent. These are used to evalu- 
ate the thermodynamical potential by adopting two different methods. First method employed the 
Berezin-Lieb inequalities to obtain an approximate form. Using this to determine the average num- 
ber of electrons (N e ) and the magnetization as well as underline their properties in terms of the the 
limiting cases (a«l,a — > oo) as well as week and strong magnetic field limits. In fact, by treating 
the limit a C 1 a correction to (N e ) is obtained. More importantly, we showed that a can be tuned 
to reproduce both the orbital paramagnetism and the Landau diamagnetism in such limit. In fact, we 
derived a general magnetic moment that can be fixed to reproduce different results and end up with 
some conclusions 

In the Second method, we employed some mathematical toy to determine the exact formula of 
the thermodynamical potential and therefore evaluated different physical quantities. More precisely, 
the Fermi-Dirac trace formulas is used and the average number of electrons as well as the magnetic 
moment are calculated. After evaluating the susceptibility, we found that there is a correction to the 
Landau diamagnetism, which is a-dependent. Again by fixing the parameter, other results can be 
obtained and in particular for a = standard results [7] is easily recovered. 

Some interesting questions remain to be solved for the present system. In fact, first concerns the 
temperature limits of the thermodynamical potential obtained in terms of the second method. This 
can also be investigated further by considering all limiting cases of the couple (a, B). Second is related 
to discuss the spatial density of current. Finally, a numerical study if the obtained results is much 
needed to give another comparisons with already published results. 

Acknowledgment 

This work was completed during A J visit to Max Planck Institute for the Physics of Complex Systems, 
Dresden. He would like to thank the Institute for the warm hospitality. He is also thankful to Mr. M. 
Said for his administrative help. 



21 



Appendix: Fermi-Dirac trace formulas 



It is well known that, like the Gaussian function, the function sechx = 1/coshx is a fixed point for 
the Fourier transform in the Schwartz space: 

I l r+oo e -ixy 

dy. (135) 



cosh s/\x sjQm 7-oo cosh \J\y 
Hence, given an Hamiltonian 7i, we can write for the corresponding Fermi operator: 

l r+oo e -(ifc+l)| (H-fJ.) 



f(H) = ^ttt^ = / — ; r^r, — dk. (136) 

JK ' \ + e P{H-v) ./_«, 4coshp v ' 

Similarly, we can write for the thermodynamical potential operator: 

— log (1 + = (2coshft)(it + 1) ^ (137) 

Therefore, the average number of fermions and the thermodynamical potential can be written (at least 
formally) as follows: 

r+oo (ife+l)& 

(N) = Tr f(H) = J ^ 4cosh|fc 9(fc)dfc (138) 

1 1 r+oo Jik+1)&£ 

n = Tr (^io g (i + e -^--)), = --y_ (2mhft)(tt + 1) ew* (139) 

where designates the function 

0(fc) = Tr(e- (ifc+1) i w ). (140) 

Observe that (2m + l)i, m G Z are (simple) poles for the function 1/ cosh Sfc and i is a pole for the 
functions @(k) and l/(ik + 1). These Fourier integrals can be evaluated by using residue theorems if 
the integrand functions <&i(k) = Q(k)/ cosh |/c and &2(k) = Q(k)/((ik+l) cosh |/s) satisfy the Jordan 
Lemma, that is, <5> x {Re i6 ) < g{R), <f> 2 (Re ie ) < h(R), for all 9 £ [0,vr], and g(R) and h(R) vanish as 
R — > 00. The quantities (iV) and O are then formally given by 



2vri 



-i(f) + a_ 1 ((2m + l)i) + ^a_ L (fc i , 



m=l 



(141) 



where a-i(-) denotes the residue of the involved integrand at pole (•), and the k u 's are the poles (with 
the exclusion of the pole i) of @(k) in the complex fc-plane. 

We now introduce the spectral resolution of the (bounded below) self-adjoint operator 7i: 

r+oo 

<p(H) = / <p(\)E{d\) (142) 



where <p is a complex-valued function and E\ = E(d\) is the resolution of the identity for the 
Hamiltonian 7i. Define the density of states u(X) as TrE(dX)/dX. The trace formula ensues: 

/+00 
tp(X)u(X)dX. (143) 
-00 
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Let us now introduce the weighted density of states w(X) = e '^ x v{\) and its Fourier transform 



-oo 



V}(k) = -^= I e~ tkA w(\) dX. (144) 



Then, from (138), (139) and (143), we can represent (N) and Q as follows: 

,+oc (ifc+l)%i Q n 

w = *L «if ?'*i tf2 '(-') < 145) 

2i /■+" e <"= +1 >* ..0,. ., 1* ft 



n - p™Eiw*+i)* ( f* ) *--r T * ,( -'' ) <146) 



where we have introduced the weighted functions 



Z x = sech(^-k)w(k), Z 2 = sech(-Trk/{3)(i2k/[3 + l) _1 u>(ife). (147) 
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